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An influence of some incorrectness of analysis on the level densities and radiative strength 

■ functions derived from the experimental 7-spectra is considered. It was shown that the ob- 

O ■ taining of reliable data from the reaction (n, 27) requires one to derive dependence of the 

^ . two-step cascade intensities on their primary transitions energy. The influence of some con- 

' ditions of an analysis of the experimental 7-spectra from reaction {^He, a) on the expected 

I value of both level density and radiative strength functions was estimated. The ways to 

C/^ [ decrease these uncertainties are suggested. 

(N 

1 Introduction 

o ■ 
a^ 
o 

Tij- ' The form of 7-spectra measured in different experiments is ambiguously determined by 

^ ' the level density p = D^^ and radiative strength functions k =< Txi > /{Ej x A'^^'^ x -Da) 
(or in equivalent form f = k x A^^^) for the 7-transitions between the levels A — > z in a 
nucleus with the mass A. The intensity of any such spectrum for a given 7-transition energy 



X 



O I is determined by folding of these two parameters of the 7-decay of a nucleus under 

pH I study. Therefore, these parameters can be determined for any nucleus only jointly from the 
solution of reversed mathematical problem. In that case, the correspondence of the shape of 
functional relation between the measured intensity of spectrum and required mean values of 
p and k to the experiment is postulated. 

In all known cases this relation includes the product p and k. That is why, any errors of, 
for example, p are completely or partially compensated in measured spectrum by correspond- 
ing errors of k. As a consequence, any method used for simultaneous determination of these 
parameters of the 7-decay by the solution of reverses problem cannot give unique functional 
dependence between p and k for given excitation and 7-transition energies. Moreover, all 
existing methods of the analysis allow one to determine with acceptable precision only the 
sum of level densities of both parities in some spin window and sum of radiative strength 
functions for dipole electric and magnetic 7-transitions. The uncertainty in level density 
partially decreases if one uses additional experimental data on the p below the excitation 
energy 1-2 MeV and spacing Dx between neutron resonances. The use of the experimental 
ratio k{El)/k{Ml) in vicinity of the neutron binding energy Bn leads to decrease in error 
of strength functions. 

Principle impossibility to determine unique quantities of p and k was directly fixed in 
method [1] used to extract these parameters from the dependence of the two-step cascade 



intensities 



following thermal neutron capture on the energy Ei of their primary transitions. This should 
be taken into account also in the analysis [2] of the primary 7-transition spectra of nuclear 
reactions. This conclusion follows from from testing [3] of possible use of the conventional 
multidimensional fitting of parameters of the Gauss method instead of the so-called "Oslo 
method" [2]. 

Multidimensional distribution of probabilities of random p and k values for the case of 
ensembles [2] of the primary 7-transition spectra, depopulating arbitrary excited level Eex, 
is degenerate even for surplus number of experimental points. As a result, equal value of 

is observed [3] for different p and k even in the case when ensemble of data [2] includes 
experimental total radiative widths, densities of low-lying levels and density [4] of neutron 
resonances extrapolated into some spin interval. 

For deformed nucleus, this extrapolation requires also the knowledge of dependence of 
coefficient of rotational enhancement of level density [4] on their spin. The lack of necessary 
experimental data on parameters of rotational bands even for the excitation energy 2-3 MeV 
and some higher results in unknown additional uncertainty in determination of p and k. Of 
course, this question is more serious for the region E^.^ ~ S„ [2]. 

Rather considerable non-linearity of the task, solution of which is described in [2] , limits 

the variation interval of p and k even at the correlation coefficient |/?corr| = 1 (as compared 
with the degenerate system of linear equations). Because the relation between the p and k for 
the primary transition spectra [2] and intensities of the two-step cascades [1] are essentially 
different, then involving of these heterogeneous experimental data in the united hkelihood 
function could increase reliability in determination of p and k. The type and magnitude 
of the most important statistics errors at both obtaining of experimental spectra in the 
frameworks of these methods and their further use for determination of p and k are partially 
separated. This also can provide increase in reliability of the results. 

Very important problem for any method for simultaneous determination of p and k is 
the minimization of additional systematic errors caused by inaccuracy in data processing. 
This follows from incompatibility of conclusions made by investigation [1] and [5] of two-step 
7-cascades following thermal neutron capture and by the study of the primary 7-transition 
spectra of nuclear reactions [6]. This is reduced to the following: according to [1], cascade 
7-decay of, for example, deformed compound nucleus is characterized by clearly expressed 
step-like structure in level density (predicted [7] by A.V. Ignatyuk and Yu.V. Sokolov and 
related by them with breaking of the next Cooper pair of nucleons). But in accordance 
with [5,6], the main peculiarity of 7-decay is the presence of "pygmy" resonance in radiative 
strength function at "smooth" level density. As the authors of [1] and authors of [5,8] use 
the same type of experimental data then all observed discrepancies should be related with 
incorrectness of a method of their analysis. This incorrectness appears to a great extent due 
to ignoring [5,8] of specific of the experiment on investigation of I^^. 



2 Specific of the experiment on investigation of the 
two-step 7-cascades 



Main properties of this experiment are the following: 

1. the impossibility to determine ordering of two successive quanta with sum energy of 
several MeV measured with HPG detectors directly in experiment; 

2. the dependence of intensity of every cascade on the p and k values over all excitation 
region of a nucleus; 

3. the connection of variation 51^^ with variations 5p and 5k through incoherent sum of 
derivatives of cascade intensities with respect to all possible values of p and k. 

1. As a result, every experimental spectrum I^^[E^) of cascades (see Fig. 1) with the 
total cascade energy Ec = Ei + E2 = Bn — Ef can be reproduce by a sum of 2^ different 
symmetrical distributions F[Ei) and S{E2) corresponding to registration of one or other 
cascade quanta as primary or secondary transitions, respectively (here are total number 
of cascades). Decomposition of the spectrum I^^{E^) into two parts corresponding to solely 
primary F{Ei) and solely secondary S{E2) = F{Ec — Ei) transitions can be realized with 
some uncertainty. Practical determination of function F{Ei) with small enough uncertainty 
is quite possible in modern experiment like [10] with the help of the nuclear spectroscopy 
methods using the numerical method [11] for improvement of energy resolution. Probable 
error of this procedure is determined by extrapolation [12] of cumulative sums of intensities 
of cascades (experimentally resolved in form of pair of peaks) for small energy intervals of 
their intermediate levels Ei < O.bBn to the zero detection threshold of their intensity. 

It should be noted that any from 2^ possible variants of functions F{Ei) and S{E2) can 
be always reproduced by some set of functional dependencies p{Egx) and k{Ei). Level density 
for any function F[Ei) is defined as the number of cascade intermediate levels in a given 
interval of excitation energy independently on the truth in determination of their primary 
transition energy. Because the value of F is always determined for any variant of quanta 
ordering in cascades then in any case are also determined and radiative strength functions. 
This means that the experimental spectra of the two-step cascades can be reproduced, in 
principle, with its experimental precision by p and k from [6] independently on the magnitude 
of systematical errors of these parameters. 

The spectra for the studied even-odd and even-even nuclei demostrate [10] relatively 
small density of peaks for the cascades with E^ < 0.5S„. This allows one to expect that the 
asymptotic value of error of function I^^[Ei) will aspire to zero as the detection threshold 
of individual cascade also aspires to zero. The portion of unobserved cascade intensity is 
determined by the detection threshold and form of distribution of deviations of the random 
7-quanta intensities from the average for any decaying level E^x- If this distribution depends 
on the structure of the level E^x then the model notions of the 7-decay must be corrected 
both in [1] and [6] . 

It is very important that the specific form of dependence of I^^{Ei) on its parameters 
allows one precisely (x^// << 1) reproduce this function with the help of p and k varying 
in narrow enough intervals of magnitudes. (Up to now, intensities of cascades F{Ei) were 



calculated [1] for several millions of pairs of different functional dependencies of p and k). 

Quite different situation is observed at direct an analysis of the experimental spectra of 
the two-step 7-cascadcs Ijj{E^) = F{Ei) + S{E2) whose amplitudes in any interval of 7- 
transition energy depend on both Ei and E2 . Any algorithm for search of random functions 
(like that used in [1]) allows one to be convinced that the precise reproduction of the sum 
of all possible values of functions F and S is possible only when the maximum values of the 
obtained p and k exceed the minimum values by a factor of several tens (figs. 2,3). It is not 
difficult to check that at the fixation of p on the ground of any additional information, the 
spectra of cascade intensities also can be exactly reproduced for all E^ by wide ensemble 
of k. Maximum value of these k is (Fig. 4, 5) 2-10 times higher than minimum values (any 
systematic errors of p shift location and distort the shape of the region of possible k values). 

It should be especially noted that at reproducing of intensities of two-step cascades in 
172 Y]-) listed in [8], "pygmy" resonance does not appear. Irremovable ambiguity of the ex- 
periment causes a necessity of transition from the analysis of the spectra I^^{E^) to analysis 
of spectra in function of the primary transitions energy I^^{Ei) only. This circumstance is 
ignored by authors of both [5] and [8] . The authors of [8] make conclusion about correspon- 
dence of the functions p^ — /{E^.^) and k'^ — (t){E^) obtained by them from 7-ray spectra of 
reaction {^He, a) to the experimental values only on the grounds that they provide repro- 
duction of the untransformed experimental distributions of two-step cascade intensities in 
limited [5,8] interval of energies of their intermediate levels. The statement about correspon- 
dence between the experimental and tested values of p and k made on the ground of equality 
jexp _ jcai ^iian ig equivalent to execution of conditions F^^p — F'^"' and S^^^ — 5"^"' for 
any energies of the cascade 7-transitions. I.e., from the multitude of possible solutions of eq. 
(2) with two unknown values of functions F and S, the authors of [5,8] choose the unique 
solution. It is wrong mathematics operation. Inequality of the experimental and calculated 
values of F and S always will mean that the p and k parameters used for calculation of cas- 
cade intensities do not correspond to their unknown experimental values. Therefore, such 
reproduction (x^// << 1) of the experimental intensities Ijj — F{Ei) +S{E2) for any ener- 
gies of cascade transitions is necessary but not enough criterion of correspondence between 
the desired and unknown experimental values of p and k. Only analysis of the experimental 
distribution I^^{Ei) with accounting for its total error allows determination of p and k with 
minimum possible systematic error. 

2. The total width of the decaying initial (or intermediate) cascade level can be presented 
as r = Fin + Tout where two items represent sum of partial widths of cascade transitions 
in limits of energy interval chosen in [5] or [8] and out of it, respectively. These sums are 
the items of functional relation (1) of I^y with p and k. If experimental distribution /^^ 
is reproduced with the use of some p and k only for the part of energies of intermediate 
levels and is not reproduced over the whole energy interval then, most probably, this partial 
correspondence should be considered as the result of mutual compensation of errors. Porter- 
Thomas fiuctuations of the primary transition widths have zero mean value and cannot lead 
to systematic discrepancy between the experimental and calculated cascade intensities. 

So, discrepancy between experimental and calculated cascade intensities observed [8] at 
low and high energies of cascade transitions are caused by the error in determination of the 
sum Tout of partial widths of primary and/or secondary cascade transitions. In turn, this 



can manifest itself only at discrepancy between the experimental and used for calculation 
values of p and k. 

3. As it is known from mathematical statistics, error of function can be presented in the 
first approach as the product of its derivative by the argument of function. For the eq. (1), 
relative errors of p and k at any energy of excitation or transition lead through incoherent 
sum of derivatives to the different relative change in I^^ simultaneously over all energy region 
of cascade transitions. If this sum is close to zero then small relative deviation of function 
J^^, for example, can be obtained only due to very wide variations of parameters or due 
to even small change in function S and necessary for this small variation of parameters of 
the calculation. In practice, in some cases [13] fixation of p makes impossible selection of k 
which provide reproduction of I^j{Ei) (especially at low energies of the primary transitions 
of cascades and large values of p). But there is no problem to achieve practically absolute 
agreement of the experimental and calculated sums F + S ior p and k with any systematical 
errors. 

3 Main sources of systematic errors in different ex- 
periments 

Discrepancy between the main properties of the p and k values obtained in [1] and [5,6] 
means the presence of systematic error in corresponding experiments. The values oi 5p and 
Sk in method [1] are determined, first of all, by the following errors of the values of the used 
function F{Ei) as: 

(a) the error in determination of absolute value (normalization of area); 

(b) distortion of the form of functional dependence I^j{Ei) on the energy of the cascade 
primary transition. 

Both errors are limited only by the experimental conditions. At the absence of these 
errors, asymptotic values of Sp and for method [1] equals about 20%. 

Infiuence of the error in determination of the absolute value of I^^{Ei) on 5p and 5k was 
studied in [13]. Most considerable result of this modeling is that the sign and magnitude 
of errors of desired parameters depend on the energy of excitation and cascade 7-transition. 
Besides, they do not depend on the error in determination of absolute values of I^j{El) 
at some these energies. This result confirms conclusion [1] that level density, for instance, 
cannot correspond to predictions of both model of non-interacting Fermi-gas [14] and model 
of constant temperature. 

In practice, the minimum possible systematic errors of p and k in any experiment and 
with the use of the data processing data [1] arc to be expected under conditions that: 

(a) the 7-ray spectrum of the thermal neutron capture was measured within the method 
[15] (or within method with equivalent precision) and capture cross-section was determined 
with minimal uncertainty; 

(b) statistics of useful events is so high that experimentally resolved peaks concentrate 



main portion of intensity of cascades with energy of intermediate levels Ei < 0.5S„; 

(c) all the other errors characterizing experiments on measurements of 7 — 7 coincidences 
were minimized. 

Because realization both of two last conditions in real experiment on measurement of 
7 — 7 coincidences is impossible then minimization of dp and 5k can be achieved at careful 
optimization of its geometry. 

Quantitative estimation of all sources of systematic errors for the data [6] was not per- 
formed. 

Systematic uncertainties of determination of p and k from the spectra h of the primary 
7-transitions from the {^He, a) reaction are mainly determined by the following: 

(a) distortion of form of the dependence of spectrum h on the primary transition energy 
at its extraction from the total experimental 7-spectra /; 

(b) degeneracy of multi-dimensional distribution of random deviations of desired param- 
eters from the expected value; 

(c) necessity of model extrapolation [2] of density of neutron resonances into the spin 
region where is no experimental data. 

These uncertainties can significantly change the conclusions on the cascade 7-decay pro- 
cess. So, Fig. 5 demonstrates estimation of interval of possible values of k{El) + k{Ml) for 
level density used for calculation shown in Fig. 3 (this density is less than that given in [6,8]). 
It maximum deviation from the previous variant is observed at Eg^ = 4 MeV and is equal to 
about 2/3 of the level density [6,8]. But even so small change in the calculated level density 
brings noticeable change in mean values of radiative strength functions which are necessary 
for reproduction of the two-step cascade intensities listed in [8]. For example, it changes 
the derivative of the strength function with respect to the energy. This occurs due to two 
factors: a necessity simultaneously to reproduce in calculation for increased level density of 
both cascade intensity and total radiative width < Fa >. The data shown in figs. 4 and 
5 point to necessity of maximum possible minimization of errors of the experiments under 
consideration. 

The authors of [2,6] do not make doubts about necessity of transition from the exper- 
imental total 7-spectra to spectra h of only primary transitions although this operation 
cannot be carried out [15] without bringing in of additional uncertainties in p and k even in 
principle. Besides, in some situations the non-linearity of the error transfer (for instance, of 
the error in determination of reaction cross-section on the error of intensity of total 7-spectra 
/) can lead even to increase in relative uncertainty of spectra h{Ei) of primary 7-transitions 
as compared with analogous value for the total spectra fi of 7-rays depopulating levels E^x- 
This conclusion follows from consideration of given in [16] and slightly modified equation 

h^fi- Y^ihifi) (2) 

i 

Here, i is the number of energy interval. The left and right parts of this equation must have 
the same dimension (events per one decay or per given number of decay in general case). 
Therefore, all spectra fi must be normalized to the same value as spectrum Fi. If it was 
done with some error, for example, fi — Ki(l)i, then it will be simply convinced that the 



re-normalization of the determined from (3) spectrum h to 100% will not have error in only 
case when Ki = const. Hence, the total spectra fi must be normalized to one decay with 
maximum possible accuracy before calculation of the primary transition spectra hi from eq. 
(3). This operation must be carried out with a precision considerably exceeding the accuracy 
in determination of the 7-ray intensities of the thermal neutron radiative capture in modern 
[15] experiment. 

Estimation of permitted Sh for the distorted spectra of any kind can be easily obtained in 
the following way: the total 7-ray spectrum / for any energy interval is recurrently calculated 
for the arbitrary p and k without statistic errors beginning from the known 7-transitions 
depopulating decay of low-lying discrete levels according to the relation 

fi^h + Y^ihJi). (3) 

i 

Then, the spectra of primary transitions arc restored from the spectra distorted by random 
error like = (1 -|- ei)fi according to eq. (3). The worst variant is the case when the error 
Cj is proportional to the number of spectrum and increases its intensity as increasing i (for 
numeration used in [16]). 

Using this technique one can obtain that maximum relative error of normalization of 
spectra to one decay must be considerably less than 1%. This follows from the fact that 
at increasing of Eg^ small value of h for low energy E.y is determined from difference of 
two much larger intensities of the total 7-spectra with very intence low energy secondary 
transitions. Moreover, the same energy of 7-quantum in spectrum / can be related with the 
sum of intensities of the first, second, third and so on cascade quanta at increasing of their 
intensities when excitation energy decreases. 

One can try to achieve appropriate accuracy using normalization of total 7-ray spectrum 
to one decay with the help of condition J2{Ej x J^) = Eg^- Such normalization excludes the 
problems [16] of different multiplicity of quanta in total 7-ray spectra for different excitation 
energy of decaying levels and for precise determination of reaction cross-section. But this 
can be done, probably, only partially. Detailed investigation of influence of energy resolution 
of detectors and functions of their response on systematic error can be required in this case, 
as well. The open question is systematic error connected with widening of spin window of 
levels excited by cascade transitions. 

Statistic error of the h data given in [2] is larger than 1%. Therefore, it is worth while to 
make extraction of the primary 7-transition spectrum for minimization of systematic errors 
by means of minimization of sums = J2{fi^^~ (^1 + Hi{hifi '^^)Y instead of mathematical 
method used in procedure [16]. 

4 Systematic errors of joint analysis of different ex- 
periments and possibilities of their decreasing 

Level densities and radiative strength functions presented in [6], in opinion of authors 
of this work, correspond to their most probable values because uncertainties of method [2] 



were completely eliminated due to fixation of level density at two excitation energies and 
the use of the experimental total radiative width of neutron s-resonances. Correctness of 
this statement depends only on the form of multi-dimensional distribution of probability of 
random values of desired parameters. 

As it was obtained in [3] , deviation distribution of the p and k values from their averages 
for the experimental data [6] is generated owing to correlation of these parameters (even at 

at the use in [6] of complementary with respect to [2] information). One can try to decrease 
this degeneracy or even eliminate it only by involving complementary information in the 
maximum likelihood function. 

Effective search for its maximum and identification of false solutions is realized now by 
means of the Gaussian method for solving systems of nonlinear equations. At any their 
degeneracy, the region of possible solutions is found within the programs of multi-parameter 
fitting with the use of the procedure of regularization of covariant singular matrix with 
obligatory variation of input values of desired parameters. The presence of exponentially 
changing parameters leads to appearance [3] of non-principle problems for the use of modern 
programs which realized the Gaussian method for determination of the most probable p and 
k values for the set of the data [6] . As the experience [3] shown, there is no grounds to wait 
principle obstacles for involving of derivatives of function I^^ from eq. (1) with respect to its 
parameters in the Jacoby matrix. 

Unfortunately, there is no any grounds to hope for complete exclusion of systematic 
errors of p and k even in this variant of analysis. Any hypotheses about the cascade 7-decay 
process of an excited nucleus (directly or indirectly used in [1], [6,8]) also lead to appearance 
of additional systematic errors of the desired parameters of the cascade 7-decay process. But 
the use of them is absolutely necessary, for instance, above the excitation energy 1-3 MeV 
of deformed nucleus due to the lack of any required experimental information in this region. 

First of all this concerns the problem of the shape of the energy dependence of the radia- 
tive strength functions of equal multipolarity and 7-quantum energy for different energies 
of decaying levels. Available data (theoretical calculation of matrix elements, experimental 
data [4] and intensities of the cascade population of levels up to the excitation energy < 1 
MeV [17] and higher [18]) show that influence of this factor on systematic error cannot be 
neglected. 

This problem common for both methods [1] and [2] . But there is and considerable differ- 
ence: in the case of the data [6], resulting value of k is a superposition of strength functions 
for 7-transitions of equal energy but depopulating levels with sufficiently different energy 
and structure of wave functions. In the case [1], difference between the energy dependences 
of strength functions of the secondary transitions and analogous dependence of the primary 
transitions affects the desired parameters p and k only indirectly - through distortion of the 
Tjf vahie. At low energy of the secondary transitions this distortion decreases due to both 
factor E'^ and possible [18] different in sign change in the shape of dependence k — <i){E^) 
for different energy of 7-transitions (for equal energy of decaying level) . 



5 Conclusion 



Insufficient degree of correctness of using mathematics and mathematics statistics for ex- 
traction of the level density and radiative strength functions from the primary 7-transition 
spectra of nuclear reactions and spectra of the two-step 7-cascades leads to mistaken conclu- 
sions about the process under study. Therefore, the use of the analysis procedure suggested 
in [19] cannot provide obtaining of the reliable data on the radiative strength functions. 
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Fig. 1. Top: the experimental spectrum of the two-step 7-cascades terminating at the first 
and ground states of "^^Al. Bottom: distribution [9] of the same intensity in function of the 
primary transition energy only. 



10% 




E , MeV 

ex' 



Fig. 2. The interval of variations of the level density providing reproduction of the experi- 
mental intensities [8] of the two-step cascades to the ground and first excited state of ^^"^Yb 
(points with bars). Line 1 represent predictions according model [14], lines 2 and 3 show the 
data used in calculation for figs. 4 and 5, respectively. 
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Fig. 3. Interval of variations of the strength functions (points 1 with bars). Points 2 and 3 
represent strength functions from [6] and [8], respectively. 




5. The same as in Fig. 3. Calculation with level density shown in Fig. 2 by line 3 



